Real-time electron dynamics with exact-exchange time-dependent density-functional 

theory 
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The exact exchange potential in time-dependent density-functional theory is defined as an orbital 
functional through the time-dependent optimized effective potential (TDOEP) method. We numeri- 
cally solve the TDOEP integral equation for the real-time nonlinear intersubband electron dynamics 
in a semiconductor quantum well with two occupied subbands. By comparison with adiabatic ap- 
proximations, it is found that memory effects in the exact exchange potential become significant 
when the electron dynamics takes place in the vicinity of intersubband resonances. 
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Time-dependent density-functional theory (TDDFT) 
P, 0] is an increasingly popular approach to the elec- 
tron dynamics in atoms, molecules and solids. Much of 
TDDFT 's success is due to its ability to describe elec- 
tronic excitations in large molecules Q using sim- 
ple time-dependent exchange-correlation (XC) function- 
al such as the adiabatic local-density approximation 
(ALDA). However, neither the ALDA nor any adia- 
batic gradient-corrected XC functional captures multi- 
ple or charge-transfer excitations 0, Q or excitons 
These difficulties are not caused by intrinsic deficiencies 
of TDDFT, but indicate the need for better function- 
al: for multiple and charge-transfer excitations one must 
abandon the adiabatic approximation, and excitons in 
solids require long-range XC kernels with a behav- 
ior. Moreover, for strong-field multiple ionization it was 
suggested that the XC potential should change discon- 
tinuously with the number of electrons N 

Several promising ideas for new time-dependent XC 
functionals have been explored, such as current-TDDFT 
[^, and many-body Green's function techniques 
Each of these approaches has its advantages and draw- 
backs: current-TDDFT is well-suited to describe polar- 
izability and collective excitations in extended systems 
[ill [ij, but introduces spurious dissipation in finite sys- 
tems|13|. The many-body techniques of Refs. give 
excellent optical spectra in insulators, but cannot be eas- 
ily extended into the nonlinear or the real-time domain. 

In static DFT, orbital-dependent functionals have long 
been recognized as powerful and versatile [l3l . If the XC 
energy is given in terms of the orbitals (e.g., the exact 
exchange energy), the associated local XC potential is 
constructed with the optimized effective potential (OEP) 
method [l5| or by the simplified but nearly as accurate 
KLI scheme [31 . This leads to XC potentials that are 
self- interaction free, have the correct — l/r asymptotics 
for finite systems, exhibit discontinuities upon change of 
N , and generally produce high-quality orbitals, eigenval- 



ues and band structures [14 1. 



This paper deals with the time-dependent optimized 
effective potential (TDOEP) [13, (Hi, which generalizes 



the static OEP method, carrying over the desirable prop- 
erties mentioned above and introducing new features that 
are unique to the dynamical case. The TDOEP integral 
equation for the local XC potential Kcca is given by [17| 
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The 0,0- are time-dependent Kohn-Sham orbitals. 
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and is an orbital-dependent XC action functional 
(which is rigorously defined on a Keldysh contour ,18]). 
Here and in the following, we formally include a spin 
index a, but consider only nonmagnetic systems. 

Linearization of the TDOEP approach leads to a 
frequency-dependent XC kernel [3], which was success- 
fully applied to optical absorption spectra of insulators 
^2Qi] and dynamic polarizabilities of atoms 2l|. In these 
applications, the frequency dependence of the linearized 
XC potential was found to play only a minor role. 

Until now there have been no applications of the full, 
real-time TDOEP scheme, only of the TDKLI approxi- 
mation [2, [l7| . A previous unsuccessful attempt at solv- 
ing Eq. ^ was plagued by numerical instabilities [13] . In 
this paper, we present an algorithm for stable numerical 
solutions of the exact-exchange TDOEP integral equa- 
tion ([T]), and apply it to the nonlinear electron dynamics 
in quantum wells. 

Our goal is to explore the significance of memory ef- 
fects in TDOEP. Memory-dependent XC potentials in 
current-TDDFT cause elastic and dissipative effects in 
the electron dynamics In particular, the adiabatic 

approximation was shown to break down in the limit of 
large, rapid deformations [23|. Here, we study the so 



far unresolved question of the importance of memory in 
exact-exchange TDDFT in different dynamical regimes. 
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To isolate the effects of memory in a clear-cut way, we 
will compare the full TDOEP with an adiabatic approx- 
imation (AOEP). 

Adiabatic TDOEP scheme. Under the implicit as- 
sumption that a system evolves so slowly that it al- 
ways remains close to the ground state of a given time- 
dependent potential, adiabatic XC functionals have no 
memory and only depend on the instantaneous density 
or orbitals. While this is straightforward in ALDA and 
TDKLI, the AOEP is more complicated since the static 
OEP 15, also depends on the energy eigenvalues, 
whose meaning in the dynamical regime is not obvious. 

We define the AOEP at time t as that static OEP 
potential whose associated ground-state density equals 
the instantaneous density n{t). This requires two steps: 
First, find that static Kohn-Sham potential V^g (e-g-j 
by iteration [23 |) which produces n{t) as its ground-state 
density, and obtain the associated complete set of orbitals 
(pj„{r) and eigenvalues e'^.. Next, take these (t)^j„{r) and 
and plug them into the static OEP integral equation 
This yields the AOEP XC potential at time t. 

TDOEP for quantum wells. We consider conduction 
electrons in rt-doped semiconductor quantum wells in 
effective-mass approximation 0, 113) El, confined along z. 
The ground-state envelope function for the jth subband 
4fj^(z) follows from a one-dimensional Kohn-Sham equa- 
tion, with density n{z) — J^'j'^'^ \'^^%i^)\'^i^F ~ <^ja)/T^ and 
subband and Fermi energy levels e^^ and ep- To describe 
intersubband dynamics preserving the translational sym- 
metry in the quantum well plane (ignoring disorder and 
phonons), we propagate the subband envelope functions 
using the time-dependent Kohn-Sham equation 
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+ ydr(2,i) +Konf(z) 



+ Vn{z,t)+Vy,^^{z,t)\^(t)j„{z,t) , (3) 

with initial condition 4'j„{z,to) = (j>'j^{z). Here, Vdr is a 
time-dependent driving field (see below), Vconi{z) and Vh 
are the confining square-well and the Hartree potential. 

Assuming that the system is in its ground state for 
t < to, one obtains the following TDOEP equation: 



= *^(fcir)' / dt' / dz'[V^„{z\t') 

j=l Jto J 
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where [k-'pY = 2(eF — eju), u^ja for quantum wells i 
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obtained following Proetto et al. |25[, and V^^, u°,^ and 



FIG. 1: (Color online) Top: dipole moment of free charge- 
density oscillations in a quantum well with two occupied sub- 
bands [Blue, red, green (dark, medium, light gray): TDOEP, 
AOEP, TDKLI]. Bottom left: TDOEP dipole power spec- 
trum. Bottom right: convergence index G (see text). 



follow from a static OEP calculation. The second 
term in Eq. ^ comes from (l)ja{z,t) = 0^g.(z)e~"^^* for 
t < to. 

Numerical algorithm. Our numerical TDOEP ap- 
proach uses a uniform spatial grid along z and a time 
discretization in uniform steps At, from to up until some 
final time T. Let us first consider two separate prob- 
lems: (i) Assuming that a V^^{z,t) is exphcitly given for 
to < t < T, Eq. ([3|) is easily propagated with the stan- 
dard Crank-Nicholson and predictor-corrector schemes 
■2]. This yields the orbitals 4>'j„{z,t). (ii) In turn, as- 
suming some (j)j^{z,tys to be given for to < t < T, Eq. 
(|4]) can be solved by discretizing the spatial and time in- 
tegrals (e.g. with the trapezoidal rule), which leads to a 
linear equation determining V^'^{z,t). 

The full TDOEP scheme requires the simultaneous so- 
lution of Eqs. (m and ^ over the interval [io,r]. We 
achieve this using a straightforward iterative loop, tak- 
ing the orbitals calculated in step (i) as input to step 
(ii), and then feeding the resulting XC potential back as 
input to step (i). Selfconsistency is reached if (p'-^ = 
and = V;'^. The loop is initiahzed with the TDKLI 
approximation for V^^ in the first iteration step. 

Let us define d = j^^dt\di{t) - di^i{t)\/ jj^dt\di{t)\ 
to monitor convergence, where di{t) — J zni{z,t)dz is 
the dipole moment at the ^th iteration. Our algorithm is 
stable and rapidly convergent, as we will show below. 

Results and discussion. We consider a 40-nm square 
GaAs/Alo.sGao.rAs quantum well with conduction band 
effective mass m* = 0.067m and charge e* = e/vT3, and 
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with an electronic density Ns ~ 2.2x lO^^cm"^ such that 
the two lowest subbands {j = 1,2) are occupied. 

We begin with the case of free charge-density oscil- 
lations. The initial state is calculated in the presence 
of a O.OlmV/nm static electric field ("tilted" quantum 
well). At to = 0, the field is abruptly switched off, which 
puts the electrons in an excited state and triggers col- 
lective charge-density oscillations. Fig. [T] shows that the 
convergence index G of our numerical algorithm drops 
down to 10~^ after only 20 iterations. Convergence was 
similar in all our TDOEP calculations. We found that 
the zero-force theorem was always satisfied to within the 
limits of numerical resolution, even for TDKLI (despite 
recent reports to the contrary in metal clusters [26|). The 
strong quantum well confinement seems to help enforce 
the zero-force theorem in TDKLI. 

Fig. [1] shows the dipole moment d{t) obtained with 
TDOEP, AOEP and TDKLI, as weU as the TDOEP 
dipole power spectrum (the others are very similar). The 
dynamics is dominated by the 1^2 and 2^3 inter- 
subband plasmons; the former dominates since 90% of 
the electrons sit in the first subband at the given Ng. 
Higher plasmons (1^4, 2^5) are orders of magnitude 
weaker. The differences between the methods are mi- 
nor: the TDOEP charge-density oscillations are slightly 
faster, TDKLI is slowest, and AOEP in between, but 
much closer to TDKLI. The first plasmon frequency 1^12 
is 10.4 meV in TDOEP and 10.2 meV in AOEP and TD- 
KLI; LU23 — 16.0 meV in all three approaches, and uju 
and W25 are at 42 and 61 meV, respectively. 

We find that the exact-exchange TDOEP does not 
cause any dissipation, similar to the high-frequency limit 
of current-TDDFT @, [l^; the nonadiabatic XC con- 
tribution is thus purely elastic, i.e. phase shifted by tt 
with respect to the adiabatic part [i^l . This is consistent 
with the observed behavior of the TDOEP versus AOEP, 
where the memory leads to a small frequency renormal- 
ization (blueshift) of the dominant UJ21 plasmon. 

Nonadiabatic effects play a pronounced role at high 
frequencies, when the system rapidly undergoes large de- 
formations. The crossover from the low- to the high- 
frequency region was found to occur around the average 
plasma frequency of the system [2^ . To explore different 
dynamic regimes, we now consider charge-density oscilla- 
tions driven by Vdi-{z,t) = e£ z f (t) sm{ojt) , with electric 
field amplitude £, intensity / ~ and frequency uj. The 
envelope f{t) is switched on at initial time Iq ~ over 
a 1-cycle linear ramp and then kept constant throughout 
the propagation. 

Figs. HOH andOshow d{t) for u = 11.2, 20, 40, and 
50 meV, probing the dynamics close to 0^12, W23 and UJ14, 
and between UJ14 and W25- The intensities are 10, 20, 40, 
and 1000 W/cm^ respectively. As expected, the dipole 
response is largest at 11.2 meV since we are close to the 
dominating plasmon. At 50 meV we need a much larger 
intensity since we are quite far away from any resonance. 
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FIG. 2: (Color online) Top: Dipole oscillations driven by an 
external field switched on at t = 0, with frequency uj — 11.2 
meV and intensity I = 10 W/cm^. Bottom: XC power [Eq. 
([5])]. Blue, red, green (dark, medium, light gray): TDOEP, 
AOEP, TDKLI. 
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FIG. 3: Same as Fig. [2 with tj = 20 meV and / = 20W/cm^. 



Again, we find that the three methods give compara- 
ble results for d{t). Figs. [2] and [3] show that AOEP falls 
in between TDOEP and TDKLI, but remains closer to 
TDKLI. For w = 40 meV there are more pronounced 
differences, probably due to the crosstalk of several plas- 
mon resonances slightly off tune. Interestingly, at w = 50 
meV, the results for d{t) are very close. We found a simi- 
lar behavior at low frequencies well below UJ12 ■ This sug- 
gests that at the intermediate frequencies we are studying 
here (not too far from the lowest intersubband plasmons), 
memory effects play a significant role only in the vicinity 
of resonances. 

For a more detailed analysis, especially of the phase 
shifts of Vxcr, it is useful to consider the power associated 
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FIG. 4: Same as Fig. E] with uj = 40 meV and / = 40W/cm^ 
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FIG. 5: Same as Fig. [1 with = 50 meV and 7 = IkW/cm^ 



ing patterns (in particular at 40 meV), which are due 
to the detuning between the driving field and the plas- 
mon resonances. These patterns look similar in AOEP 
and TDKLI, but exhibit marked additional structures in 
TDOEP. Remarkably, as seen most clearly at 50 meV, 
these strong and rapid fluctuations of V^a in TDOEP 
leave hardly any imprint on d{t). 

Conclusion. We have analyzed the electron dynam- 
ics in quantum wells with exact-exchange TDOEP, us- 
ing a stable and rapidly convergent numerical scheme. 
Away from resonances, the AOEP and TDKLI closely 
agree with each other (like ground-state OEP and KLI 
[Iq ) and are good approximations to the full TDOEP. 
Memory effects become more significant in the vicinity 
of intersubband resonances, resulting in additional elastic 
contributions to the dynamics. Exact-exchange TDOEP 
has no memory at all in systems with only one occu- 
pied level [16|], and our quantum well has only a small 
population of the second subband. We expect nonadia- 
batic effects to become more important in systems with 
a larger relative occupancy of upper levels, and in fre- 
quency regimes extending further beyond the lowest ex- 
citations. 

Our analysis of the dynamic XC force shows that the 
exact-exchange V^a has much richer temporal features 
than the adiabatic approximations, but this has relatively 
little impact on the electron dynamics itself. This sug- 
gests that for intrinsically nonadiabatic effects such as 
multiple excitations ^] and dissipation [lo| one needs to 
go beyond exact exchange. 

This work was supported by Research Corporation and 
by NSF Grant No. DMR-0553485. We thank Paul de 
Boeij for valuable comments. 



with the dynamic XC force: 

P{t) = I dzj{z,t)V, [V.a{z.t) - Vl{z)\ , (5) 

where j{z,t) = YZ" '^WkA^^'^)'^ ^'t^k^{^^'^)WF - <^ka)h 
is the current density. Since there is no dissipation, P{t) 
is zero on average, but fluctuates more rapidly than the 
charge-density oscillations, at least twice as fast as d{t). 

The bottom panels of Figs. HlHlshow P{t) for w = 11.2, 
20, 40, and 50 meV. The differences between TDOEP 
versus AOEP and TDKLI are now much more appar- 
ent, in particular for higher frequencies. At 20 meV one 
sees that I^o- picks up a signiflcant phase shift compared 
to AOEP and TDKLI, which themselves are completely 
in sync. This is a clear indication of the elasticity in- 
duced by the memory. The effect becomes even more 
pronounced at 40 and 50 meV. 

We also find that P{t) exhibits characteristic beat- 
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